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Abstract 

We examine the band spectrum, and associated Floquet-Bloch eigensolutions, arising in 
a class of three-phase periodic checkerboards. On a periodic cell [—1, 1[^, the refractive in- 
dex is defined by ra^ = 1 + gi{xi) + g2{x2) with gi{xi) = for < Xj < 1, and gi{xi) = 
for — 1 < < where is constant. We find that for > — 1 the lowest frequency 
branch goes through origin with linear behaviour, which leads to effective properties en- 
countered in most periodic structures. However, the case whereby = —1 is very unusual, 
as the frequency A behaves like \/k near the origin, where k is the wavenumber. Finally, 
when < —1, the lowest branch does not pass through the origin and a zero- frequency 
band gap opens up. In the last two cases, effective medium theory breaks down even in the 
quasi-static limit, while the high-frequency homogenization [Craster et al., Proc. Roy. Soc. 
Lond. A 466, 2341-2362, 2010] neatly captures the detailed features of band diagrams. 
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Resume 

Sur une classe d'echiquiers a trois phases aux proprietes effectives singulieres 

Nous etudions le spectre de bande associe aux modes de Floquet-Bloch dans une classe 
d'echiquiers periodiques. Sur une cellule de base [— 1,1[^, I'indice de refraction est defini 
parn^ = 1 + gi{xi) + g2{x2) ou gi{xi) = (uneconstante), pour < Xj < 1, et gi{xi) = 
pour — 1 < Xj < 0. Pour > —1, la premiere bande passe par I'origine avec un 
comportement lineaire, ce qui conduit a des propriete effectives rencontrees dans la plupart 
des structures periodique. En revanche, le cas = —1 est moins ordinaire, puisque la 
bande de frequences acoustiques A se comporte comme au voisinage de I'origine, avec k 
le nombre d'onde. Finallement, quand < —1, la bande acoustique disparait: la premiere 
bande ne passe plus par I'origine et une bande interdite a frequence nulle apparait. Dans 
ces deux derniers cas de figure, la theorie des milieux effectifs ne s'applique pas, alors que 
la theorie d'homogeneisation hautes frequences [Craster et al., Proc. Roy. Soc. Lond. A 
466, 2341-2362, 2010] reproduit avec precision les diagrammes de bandes. 
Mots- des : Homogenisation ; Refraction negative ; Bande acoustique 
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Version frangaise abregee 



Nous considerons un echiquier periodique dont la geometrie est decrite dans la figure 
1. Les ondes optiques ou acoustiques qui s'y propagent sont solutions de I'equation de 
Helmholtz (1) dans les cases homogenes occupees par des materiaux isotropes qui sont 
positifs, nuls, ou negatifs. Par ailleurs, des conditions classiques de continuite de la com- 
posante tangentielle du champ electromagnetique (ou du champ de deplacement et de la 
composante normale du stress) sont imposees aux interfaces entre les differents materiaux 
de la cellule elementaire et des conditions de Floquet-Bloch sur les bords opposes de celle- 
ci. Dans le cadre de la theorie d'homogeneisation hautes frequences developpee dans [1], 
nous nous interessons plus particulierement au cas oil les conditions de Floquet-Bloch 
se reduisent a des conditions periodiques ou anti-periodiques, ce qui permet de deduire 
I'equation homogeneisee (4) qui conduit a une estimation fine des courbes de dispersion a 
travers la frequence effective (5), pour des nombres d'onde au voisinage des points A, B et 
C de la zone de Brillouin, cf. figure 1. II s'avere que I'on pent meme deduire les relations 
de dispersion analytiquement dans le cas d'echiquiers a trois phases, donnees par (7) et 
(6). 

Grace a ces approches de type homogeneisation hautes frequences [1] , et Kronig- Penney 
[5], nous mettons en exergue trois types de comportements effectifs pour des ondes se 
propageant dans de tels echiquiers aux basses frequences: ordinaire de type cristal pho- 
tonique (r^ > —1: bande acoustique lineaire a I'origine, voir figure 3) (a)), singulier de 
type cristal plasmonique (r^ = —1: bande acoustique non-lineaire a I'origine, voir figure 
2(a)) et de type metamateriaux (r^ < —1: perte de bande acoustique i.e. bande interdite 
a frequence nuUe, voir figure 2(b)). 

Nous concluons notre etude par une application de type lentille echiquier par refraction 
negative quand > —1, avec I'image d'un point source a la frequence donnee par 
I'estimation (8), voir figure 3(b). 



1. Introduction: Beyond effective medium theory with high-frequency homogenization 

Negative refraction is an emerging field in photonics introduced by Victor Veselago in 
the late 1960's [2], with renewed interest following the controversial claim of John Pendry 
that negative refraction makes a flat lens with unlimited resolution possible [3]. The quest 
for this superlens has fuelled research in structured photonic materials during the last 
decade. These composites are known as photonic crystals when the wavelength of light is 
of the same order as the typical heterogeneity size (Bragg frequency regime), and metama- 
terials when the wavelength is much larger (low frequency homogenization regime). The 
propagation of waves with anti-parallel group and phase velocities, an essential ingredient 
for negative refraction, was discussed in the late 1940's by Leon Brillouin [4]. 

In order to investigate the properties of photonic crystals, we consider the Helmholtz 
equation (1) in a periodic medium. This equation could, with appropriate notational and 
linguistic changes, hold for acoustic, electromagnetic, water or out-of-plane elastic waves 
and so encompasses many possible physical applications. We solve 
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+ + gi{xi) + g2{x2)]u = 



(1) 



+ 



(9X2^ 



for u{xi,X2) on the square —1 < xi,X2 < 1, where is the frequency squared. In the 
case of a three-phase checkerboard with a square cell as shown in Fig. 1, gi{xi) is taken 
to be the piecewise constant 



We note that in the context of optics, the unknown u in (1) stands for the longitudinal 
component of the electric field and is associated with bj^j(? whereby to is the 
electromagnetic wave frequency and c is the speed of light in vacuum. Moreover, Sry^Q = 
1 + gi{xi) + (72(2^2) where fio = I and Er are the relative permeability and permittivity of 
the dielectric (non-magnetic) medium. 

For waves through an infinite, perfect, doubly periodic checkerboard one can invoke 
Bloch's theorem [4] and simply consider the square cell with quasi-periodic Bloch boundary 
conditions applied to the edges: 

m(1,X2) = e"'iM(-l,X2),Mxi (1,3:2) = e"'iM^i(-l,X2), 
u{xi, 1) = e"^^u{xi, -1), Ux^{xi, 1) = e'^'^u^ciixi, -1), 

that involves the Bloch wave- vector k = {ni, k,2) characterizing the phase-shift as one 
moves from one cell to the next. There is also continuity of u,Uxj^,Ux2 along xi = and 
X2 = 0. This Bloch problem is solved explicitly and dispersion relations that link the 
frequency and Bloch wavenumber are deduced; as is well-known in solid state physics 
[4] only a limited range of wavenumbers need be considered, namely the wavenumbers 
along the triangle shown in the reciprocal Brillouin lattice of Fig. 1 and these form the 
irreducible Brillouin zone. The dispersion curves reported in Fig. 2-3 illustrate several 
interesting features: stop-bands for which wave propagation is not possible, and regions 
of flat dispersion curves for which the group velocity is zero and features of slow sound 
or light occur. In this paper, we focus our attention to the lowest curves, associated with 
averaged properties of the checkerboards, and discuss the extension of existing effective 
theories to higher frequencies. 

Conventional homogenization is widely assumed to be ineffective for modelling photonic 
crystals as it is limited to low frequencies when the wavelength is long relative to the 
microstructural lengthscales. Here the recently developed high frequency homogenization 
theory [1], which is free of the conventional limitations, is used to generate effective partial 
differential equations on a macroscale, that have the microscale embedded within them 
through averaged quantities. We focus our attention on periodic checkerboards, some of 
them with sign-shifting coefficients [13,15]: the latter are known to support a host of sur- 
face plasmons associated with highly resonant features less than conducive for any effective 
theory. Our aim is actually to push the newly introduced high-frequency homogenization 
beyond its limits. 

At leading order, the field m(x) propagating inside the periodic checkerboard can be 
represented as U{^)f(X.) where ^ and X are the short- and long-scales [1]. The function 
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1 < Xi < . 
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Figure 1. Left: A single cell for the checkerboard geometry, for piecewise constant nredia. Right: The reciprocal Brillouin lattice in 
wavenumber, k, — {ki^k^). space 

U represents the local short-scale solution at a specific standing wave frequency, Aq (which 
is in the high frequency regime), and / the long-scale variation. Ultimately an homogenized 
PDE emerges for / entirely on the macroscale [1] 



T- 

■J- 7,1 



dX,dX, 



+ (A^ - \l)f = 0, 



(4) 



where the spatially constant tensor Tij incorporates the short-scale information associated 
with the standing wave frequency Aq; A is the full frequency. The PDE can be augmented 
by additional terms if there is material variation leading to localised defect modes [7] . The 
PDE (4) can then be utilized for, say, Bloch waves where / ~ exp(ik ■ X/2), to find local 
dispersion relations 



A ~ A 







2 »J j 



(5) 



valid near the standing wave frequency Aq; this formula is explicitly for those periodic- 
periodic standing waves at wavenumber A and almost identical formulae hold for the other 
standing waves. Thus, the constant T^j, once identified, completely encapsulate the effect 
of the microstructure on the dispersion properties of the system; the disparity of scales is 
encapsulated in e which is assumed positive and <C 1. Some typical values of the tensor 
Tij are given for the cases = —1 and —2 in table 2. One should note that in some cases 
the diagonal entries, are identical (effective isotropic parameters), different (anisotropic 
parameters) and even of opposite signs (a hallmark of anomalous dispersion). 



2. An exact dispersion relation for a class of three-phase checkerboards 

We have first investigated the dispersion curves associated with three-phase periodic 
checkerboards using the finite element software COMSOL MULTIPHYSICS. However, 
the convergence of the numerical algorithm can be hazardous, due to sign-shifting coef- 
ficients [15]. Thankfully, the Bloch problem for such checkerboards can be solved ana- 
lytically and this provides probably the only non-trivial two dimensional structure with 
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= -1 = -0.9 

A /i Til T22 A li Til T22 

AO N/A N/A AO JV/A N/A 

A 4.801 2.605 -2.907 0.195 A 4.678 2.544 -4.436 0.278 

A 4.801 2.605i 0.195 -2.907 A 4.678 2.544i 0.278 -4.436 

C 2.652 -1.813 -1.813 C 2.6086 -2.1778 -2.1778 

C 4.519 2.345 2.951 -0.380 C 4.263 2.1596 4.079 -0.678 

C 4.519 2. 345i -0.380 2.951 C 4.263 2.1596i -0.678 4.079 

B 2.336 0.767i -3.445 1.394 B 2.2473 0.8165i -4.3226 1.4607 

B 4.477 2.350i 3.215 0.384 B 4.1814 2. 1623i 4.7095 0.6936 

B 4.821 2.595 -0.199 -2.774 B 4.708 2.530 -0.2883 -4.1548 



Table 1 

The standing wave frequencies, A, the /i and the corresponding Tn and T22 for tlic lowest three dispersion curves for — —1 and 
= -0.9 



A ^l, Til T22 

A 2.8898 0.4519 0.4519 
C 3.0217 -0.3412 -0.3412 
B 2.9632 0.3524 -0.41 0.3596 



Table 2 

Values of A and the corresponding /i at the wavcnunibcrs for standing waves with the associated values of T\\ and T22 for the lowest 
dispersion curve for — —2. 

an exact solution; it is a natural generalization of the classical Kronig-Penney [5] one- 
dimensional piecewise constant periodic medium. The dispersion relation follows from 
setting m(xi,X2) = Qi{xi)Q2{x2) in the governing equation (1). From the piecewise con- 
stant fonction (2), and from the Bloch conditions (3), the following coupled equations 
emerge [6] 

27j/3j(cos Ki - cos 7j cos A) + {(3^ + jf) sin 7^ sin f3i = 0, (6) 

for i = 1,2. These are dispersion relations for A^,/i^ (the latter is a separation constant) 
in terms of ki, 1^2] the parameters 7^, /3j act to couple the equations as 

A2 = AV2T/i^ 7.^ = A2(l/2 + r2)T^^ (7) 

for i = 1,2 with the minus, plus signs for i = 1,2 respectively. This two variable disper- 
sion relation is solved via matrix Newton iteration. For the exact solution found in this 
section, one finds A, fi explicitly and some sample eigenvalues are given in table 2. Notably 
some frequencies. A, have two associated fi values, one real and one imaginary and these 
correspond to degenerate cases where the dispersion curves touch. 



3. Illustrative electromagnetic paradigms whereby effective medium theory fails 

Effective properties of periodic checkerboards are of particular mathematical interest, 
and their study goes back a long way, with the seminal paper by Keller [8] for two- 
phase checkerboards, and those of Mortola-Steffe [9], Cr aster- Obnosov [10] and Milton 
[11] for four-phase checkerboards. Physicists have further analysed the effective proper- 
ties of two-phase checkerboards with continuously varying parameters in the context of 
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Figure 2. Dispersion diagrams for — —1 (a) and — —2 (b). Panel (a) shows the lowest three exact dispersion curves from the 
explicit dispersion relations and the asymptotic dispersion curves with the behaviour at the lowest curve that touches the origin 
shown dotted. The lower inset shows the detail near the origin comparing the exact (dashed) and asymptotic (dotted) solutions. 
The upper inset shows a similar comparison (asymptotics given by the solid line) for the lowest curve near wavenumber position C. 
Panel (b) shows that there is no acoustic band and that the dispersion curves are nearly flat. The inset shows an enhanced view of 
the lowest dispersion branch with the numerics (dashed) shown together with asymptotics (solid). 



photonic crystals [12]. In this note, we investigate the case of three-phase checkerboards 
with negative and vanishing refractive index (when takes negative values). Such a 
study is motivated by the fabrication of nano-scale gold and silver checkerboards whose 
cells alternate positive and negative refractive index media in the visible range of fre- 
quencies. Potential applications lie in extra-ordinary transmission of light through the 
sub- wavelength aperture holes in such checkerboards, and tremendously enhanced local 
density of states for light confinement [13] . These physical phenomena are underpinned by 
the effective properties of such checkerboards. Hence, an adequate homogenization model 
for such structures is of pressing importance. 

We start to explore the case of a checkerboard with = —1. Dispersion diagrams re- 
ported in figure 2(a) compare the solution from the exact dispersion relations with those 
from the asymptotics. It clearly demonstrates the superiority of the high-frequency ho- 
mogenization approach which reproduces very precisely the acoustic and optical branches, 
thereby extending classical homogenization to the stop band (Bragg) regime of frequen- 
cies. Moreover, the region where the lowest branch cuts the origin, the long wave low 
frequency regime, is no longer linear and asymptotics of the dispersion relation show that 
the local behaviour is that A ~ 6^/^|k|^/^. The insets show detail of the asymptotics versus 
the exact solution. 

Figure 2(b) shows the case of a checkerboard with = —2; the acoustic band is lost 
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for < —1, which is reminiscent of singular problems in homogenization of periodic 
arrays of infinitely conducting fibres in transverse electric polarization (modelled with 
Dirichlet boundary conditions) [14]. In the present case, we do not consider Dirichlet 
conditions (corresponding to hard walls in acoustics) but we rather have sign-shifting 
conditions across some interfaces. It is clear that effective medium theory breaks down, 
however the high-frequency homogenization approach captures the fine features of the 
nearly flat dispersion curves, see the inset to figure 2(b), which are associated with slow 
waves, another topical subject in photonics. 

We finally study the case of a checkerboard with = —0.9, figure 3(a) looks visually 
similar to figure 2(a); however, a closer view on a log- log scale (not shown) of the acoustic 
band reveals that the long wave low frequency limit differs from that when = — 1. 
For = —0.9 one gets the usual effective medium result with frequency linearly related 
to wavenumber however for = — 1 the effective medium theory no longer holds and 
a different asymptotic relation holds (A ~ 6^/^|fi;|^/^); we further observe that the high- 
frequency homogenization approach (dotted line) captures the fine features of this curve 
(dashed line). 

Figure 3(a) also shows the light-line emerging from C relevant for waves impinging upon 
a checkerboard structure rotated through 7r/4. Notably the light-line crosses the lowest 
dispersion curve in such a way that their group velocities are opposite and therefore one 
can induce negative refraction at this frequency. The asymptotic theory explicitly identifies 
this frequency as 



A = Ao/^l - T/A (8) 

where Tn = T22 = T for that case. The estimate is that A ~ 2.1 and computations for 
this are shown in figure 3. 




Figure 3. Dispersion diagrams for — —0.9. Left panel shows the dispersion curves that are visually close to those for — — 1 
however, the behaviour near A for low frequency is asymptotically different between these cases. The dot-dash line is the light-line for 
waves impinging upon the rotated checkerboard. The right panel shows the response due to a harmonic point source at frequency 
A — 2.1 above a finite size checkerboard for — —0.9 with each cell rotated by an angle 7r/4; it displays an image below the 
checkerboard in accordance with the Snell-Descartes law for a negative refractive index. 
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4. Conclusions 



We have investigated the Bloch spectrum of periodic, three-phase, checkerboards. The 
focus has been on special cases in which the refractive index squared (i.e. the relative 
permittivity in the contex of optics) in homogeneous cells can take vanishing and/or neg- 
ative values, and we have found that anomalous dispersion effects of waves are possible, 
including acoustic bands with non linear features for vanishing wavenumbers, zero fre- 
quency stop bands, and nearly flat dispersion curves. The last two items are known to 
be respectively associated with singularly perturbed problems in electromagnetism that 
cannot be homogenized in a conventional way [14], as no wave is allowed to propagate 
in the periodic structure at low frequencies, and slow waves, whereby delay lines can be 
achieved. 
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